Mott metal-insulator transition in the Hubbard model 
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The ground state of the Hubbard model is studied within the single-site approximation (SSA) 
and beyond the SSA. Within the SSA, the ground state is a typical Mott insulator at the critical 
point n = 1 and U/W = +oo, with n being the electron density per unit cell, W the bandwidth 
of electrons, and U the on-site repulsion, and is a normal Fermi liquid except for the critical point. 
Beyond the SSA, the normal Fermi liquid is unstable against a non-normal Fermi liquid state 
except for a trivial case of U — such as a magnetic or superconducting state in two and higher 
dimensions. In order to explain actual observed metal-insulator transitions, one or several effects 
among the electron-phonon interaction, multi-band or multi-orbital effects, and effects of disorder 
should be considered beyond the Hubbard model. 
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The Mott metal-insulator (M-I) transition is an inter- 
esting and important issue in solid-state physics [lj , and 
a lot of effort has been made towards clarifying it Q. 
However, its theoretical treatment is still controversial. 
One of the most contentious issues is whether or not the 
transition can be explained within the Hubbard model. 

In Hubbard's approximation 0, 0], a band splits into 
two subbands when the on-site repulsion U is large 
enough such that U > W, with W the bandwidth 
(W > 0); the subbands are called the upper Hubbard 
band (UHB) and the lower Hubbard band (LHB). In 
Gutzwiller's approximation 5, 6], a narrow quasiparticle 
band appears around the chemical potential; the band 
and quasiparticles are called the Gutzwillcr band and 
quasiparticles. It is plausible to speculate that the den- 
sity of states in fact has a three-peak structure, with the 
Gutzwillcr band between UHB and LHB. Both of the ap- 
proximations are single-site approximations (SSA). An- 
other SSA theory [7[ confirms this speculation, showing 
that the Gutzwiller band appears at the top of LHB when 
the electron density per unit cell n is less than one, i.e., 
n < 1 . According to Kondo-lattice theory 0, 0, [l(| , the 
three-peak structure corresponds to the Kondo peak be- 
tween two subpeaks in the Anderson model, which is an 
effective Hamiltonian for the Kondo problem. An insu- 
lating state can only be realized when the Fermi surface 
(FS) of the Gutzwiller quasiparticles vanishes. 

When n = 1 and W/U = 0, a localized electron or a 
non-interacting spin sits at every unit cell. In this case, 
the ground state is infinitely degenerate and is a typi- 
cal Mott insulator. In Gutzwiller's approximation, the 
ground state is a metal for n =^ 1 or W/U > 0. Accord- 
ing to Brinkman and Rice's theory [ill ] . when n = 1 the 
effective mass of the quasiparticles diverges at a critical 
U ~W, which is denoted as Ubr here. It may be argued 
that the FS of the quasiparticles vanishes for n — 1 and 
< W/U < W/Ubr, i.e., on the dashed line in the phase 
diagram shown in Fig. [TJ It is curious that no order pa- 
rameter appears in this second-order transition and every 
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FIG. 1: Schematic phase diagram for T = K within the sub- 
space with no order parameter. The arc indicates a possible 
but unlikely first-order transition line between metallic states. 
Dots indicate critical points. See the text for the dashed line. 



physical property seems to be continuous across the line. 
It is probable therefore that the critical U is infinite be- 
yond Gutzwiller's approximation. One of the purposes 
of this Letter is to show that no Mott M-I transition is 
possible at any finite U, which contradicts a scenario that 
actual observed M-I transitions can be explained within 
the Hubbard model. Hence, the other purpose is to ex- 
amine relevant effects beyond the Hubbard model, which 
should be considered to explain the transitions. 
The Hubbard and Anderson models are defined by 
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unit cells. When it is assumed there is no order parame- 
ter, the self-energy of electrons for the Hubbard model is 
divided into single-site and multi-site ones: E CT (ie n ,k) = 
S<j(i£ rl ) + AS (T (ie„, k). The single-site one S ff (i£„) is 
given by that of the Anderson model when U and the 
single-site electron lines of the Feynman diagrams are the 
same in the two models. The chemical potentials of the 
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two models are denoted by p and p, respectively. The 
single-site Green function of the Hubbard model is given 
by R a (ie n ) = (l/JV)X)k 1 /[en + M--B(k)-S (r (ie n) k)], 
with E(k) = e„ - (l/iV)X; <#j ^-exp[ik- (Ri-Rj)j. 
Here, Ri is the position of the ith lattice site. The Green 
function of the Anderson model is given by G a {ie n ) — 
l/[ie„ + p - e<j - Soften) - L CT (ie„)], with L a (ie n ) = 
0-/ N ) Ek I^( k )l7 i i£ n + P- E c (k)}. The condition for 
the electron lines is R a (ie n ) = G a (ie n ). In fact, a set of 
A _ e d = M - £<j and 

L ct (e-H0) = £ + /i - e Q - £ CT (£ + iO) - (e+zO) , (3) 

is a mapping condition to the Anderson model. First, the 
multi-site self-energy is ignored in Eq. ([3]). The approx- 
imation is the best SSA, which considers all the single- 
site terms, and is reduced to a problem of determining 
and solving the Anderson model [|, @, G3- The SSA 
is rigorous for infinite dimensions within the Hilbert sub- 
space with no order parameter [12j . The SSA can also be 
formulated as the dynamical mean-field theory (DMFT) 
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mation 

According to Yosida's perturbation theory [l6j and 
Wilson's renormalization-group theory (l7| . provided 
that the FS of conduction electrons is present, the ground 
state of the s-d model is a singlet or a normal Fermi liq- 
uid (FL). Exceptionally, it is a doublet for J s -d = 0, with 
J s -d the s-d exchange interaction. Since the s-d model 
is derived from the Anderson model, the ground state of 
the Anderson model is a singlet or a normal FL, or ex- 
ceptionally a doublet for the just half filling and infinite 
U. A necessary condition for the presence of the FS is 



Im [L a (e + i0)] e 



< 0. 



(4) 



This condition is called the FS condition in this Letter. 

The density of states for U = is given by po{e) = 
(l/A^)Y^ k (5[e— E(k)]. It is assumed that the FS is present 
for U = or that < n < 2 and there is no gap in pg (e) . 

As a preliminary, a Lorentzian model is examined, 
where p (e) is given by p L (e) = (l/n)A/[(e - e a ) 2 + A 2 }, 
with A > 0. According to Eq. L a (e + i0) = -iA 
[13] ]. Since the FS condition ([4]) is satisfied for the An- 
derson model, the ground state of the Hubbard model is 
a normal FL except for n = 1 and W/U = 0. 

In order to examine a non-Lorentzian model of p n {£), 
the following model is first examined: 
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with 5 > 0. According to Eq. ([3]) and the inequality 
J de' p {e')[Si{e,e') + x] 2 /[Sle,s') + S 2 (s)} > for any 
real x, with Si(s,e') and £2 (e) being real functions de- 
fined by Si(e, e') + 182(e) = £ + p — e a — £ ; — S CT (£ + iQ), 
it follows that lmL a (e + iO) < 0. The equality is only 



possible provided that S = +0. The FS condition (HI) 
is satisfied for 5 > 0. The ground state for p n (e) is 
obtained in the adiabatic process S — ► +0. This gives 
a normal FL except for n = 1 and W/U = 0, which 
may be degenerate with a non-normal FL provided that 
Im[L CT (£ + z0)] e _ > ±o.(5->+o = 0. No symmetry breaking 
occurs in the adiabatic process. 

Since symmetry breaking is caused by Weiss mean 
fields, which are intersite effects, symmetry breaking 
is definitely impossible in the SSA, as is shown above. 
Under the SSA, therefore, the adiabatic continuation 
[HI as a function of U holds. The ground state for 
the non-Lorentzian p n {e) is the normal FL obtained in 
the adiabatic process. The self-energy is expanded as 

S CT ^0) = E CT (O) + (l-0 7 )£+(l-0 s )^CT 5 MB-ff + O(£ 2 ) 

at T = K in the presence of an infinitesimally small Zee- 
man energy gpsH, with E CT (0), 6^, and <f> s all being real. 
According to the FS sum rule [19( the electron density 
n is given by n - (1/N) £ k(T 9([p- E(k) ~ t a {Q)]/W) , 
with 6{x < 0) = and 6{x > 0) = 1. The density of 
states is given by p(e) = —(l/ir)lmR a (e + iO). When 
n is kept constant, p — S ff (0), p(0), and L a (+i0) do 
not depend on U. Since ImL a (+i0) < for U = 0, 
Imi CT (+i0) < for U/W > 0. When p n {e) is continuous 
and finite, L a {e-\-iQ) is continuous and the FS condition 
is satisfied so that the FL obtained in the adiabatic pro- 
cess is not degenerate with a non-normal FL. Provided 
that is discontinuous or divergent, degeneracy is 

possible for certain n corresponding to the discontinuity 
or divergence of p n (s) according to the FS sum rule. 

In Gutzwiller's approximation, 7 oc 1/|1 — n\ for 
W/U — 0. Since <^ 7 diverges as n — > 1, Re[E cr (£-(-z0)] rl _ > i 
is at least discontinuous at e — so that Imp^e + 
i0)] e _»±o.n^i diverges logarithmically according to the 
Kramers-Kronig relation. Then, p(e — > ±0) = for 
n -> 1 and W/U = 0, although p(0) > for n = 1 
and W/U — * 0. There is a discontinuity in p(e) as a 
function of e, n, and W/U . The critical point of n = 1 
and W/U — is unconventional within the SSA. 

If the mapping to the Anderson model is not unique 
or physical quantities such as 7 and <p a are multi-valued 
functions of U, a first-order transition between metal- 
lic states occurs. The adiabatic continuation still holds, 
for example, along a route around one of the critical 
points at the ends of the first-order transition line. Since 
the FS sum rule, p(0), and L a (+i0) are all the same in 
the two metallic states on different sides of the line, the 
occurrence of the first-order transition is unlikely. The 
transition never occurs in the Lorentzian model because 
the mapping is unique. The transition is shown on a 
schematic phase diagram in Fig. [T] 

The occurrence of a first-order M-I transition is sug- 
gested by numerical theories using the Monte Carlo 
method ^ and DMFT [2l|. In these theories, the com- 
pressibility Xc(0) = dn(p)/dp shows a rapid change 22 1. 
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When the rapid change is really a jump, the phase dia- 



gram for T > K is like that shown in Fig. 1 of Ref. [21 



The phase diagram suggests that the first-order M-I tran- 
sition occurs even at T = K. If this is the case within 
the Hilbert subspace with no order parameter, however, 
it cannot be an M-I transition, but rather it must be a 
transition between metallic states. 

The irreducible spin polarization function is also di- 
vided into single-site and multi-site ones: 7r s ('w;,q) = 
7T s (iu)i) + An s (iu}i,q). The single-site one 7r s ('w/) is 
given by that of the Anderson model. The spin 
susceptibilities of the Anderson and Hubbard models 
are given, respectively, by x s (iuJi) = 2jr s (iuJi)/[l - 
Uws(ivi)] and x,(io;j,q) = 27r a (iu;j,q)/[l - Uir s (iu>i, q)]. 
In Kondo-lattice theory, an intersite exchange interac- 
tion I s (iuji,q) is defined by Xs(*^i,q) = Xs(ivi)/[1 - 
(l/4)J,(iw,,q)x.(iwi)]. When U/W > 1, I a (iw,,q) ~ 
2U Att s (iuji , q) , where 0[1/Ux s (iui)] terms are ignored. 
The exchange interaction is composed of three terms 23 1 . 
The first is the superexchange interaction. According to 
field theory, it arises from the exchange of a pair exci- 
tation of electrons between LHB and UHB. When it is 
assumed the widths of LHB and UHB are vanishingly 
small, the strength of the interaction between nearest 
neighbors is J = — 4\t\ 2 /U, with t the transfer integral be- 
tween nearest neighbors [24| ■ The second is an exchange 
interaction arising from the exchange of a pair excitation 
of the quasiparticles. The strength of the interaction is 
proportional to the bandwidth of the quasiparticles. It is 
antiferromagnetic when the nesting of the FS is sharp or 
the chemical potential lies around the band center of the 
quasiparticles, and it is ferromagnetic when the chem- 
ical potential lies around one of the band edges of the 
quasiparticles. The third term corresponds to the mode- 
mode coupling term in the self-consistent renormaliza- 
tion theory of spin fluctuations 25( , which is relevant for 
U/W < 1. According to the Ward relation 26], the static 
component of the single-site irreducible vertex function in 
spin channels is given by X s — 4> s [l — Un s (0)]. When this 
is used, the mutual interaction between the quasiparticles 
is given by {U\ s ) 2 [x s (iu>u q) - Xs(m)] = ^I*(m,q), 
with I*(iuji,q) = 7 s (iw;,q)/[l - (l/4)I s (iuJi,q)xs(m)]; 
the single-site term is subtracted because it is considered 
in the SSA. The mutual interaction mediated by spin 
fluctuations is essentially the same as that due to the 
exchange interaction I s (icui,q) or J* (icjj , q) . In Kondo- 
lattice theory, an unperturbed state is constructed in the 
non-perturbative SSA theory [27| and intersite effects are 
perturbatively considered in terms of J s (icJ;, q) [28j ]. 

When the Fock-type self-energy due to the superex- 
change interaction J is considered, the bandwidth of the 
quasiparticles is given by W* oc (W/<py) + c.j\J\, with 
c,j = 0(1) 0. When this renormalization is considered, 
the critical point of n — 1 and W/U — turns conven- 
tional; there is no discontinuity in p(e) as a function of 
s, n, and W/U. 



An order parameter can appear in two dimensions or 
higher. The stability of a normal FL against an ordered 
state with an order parameter can be examined when the 
corresponding response function is perturbatively consid- 
ered in terms of I s (iu;;,q) or I*(iuji,q). When J s (iw;,q) 
is strong, for example, the FL is unstable against a mag- 
netic state. The Neel temperature TV is defined as the 
highest value of T/v determined by [xs(0, q)]T=T N ~~ * +°° 
as a function of q. When J s (icJ;,q) is so weak that 
[Xs(0> q)]t=ok < +oo for any q, the FL is stable against 
any magnetic state. 

The energy scale of local quantum spin fluctuations 
is defined by &b7k — 1/[x s (0)]t=ok, where Tk is the 
so-called Kondo temperature. In accordance with the In- 
dependent crossover between a localized spin for T ^> Tk 
and a normal FL for T <C Tk in the Kondo problem 
fl7j |. magnetism for Tn 3> Tk is characterized as local- 
moment magnetism and magnetism for Tn <C Tk is char- 
acterized as itinerant-electron magnetism 29] . The mag- 
netism crossover is simply a Mott M-I crossover. 

When I*(iivi,q) is weak or strong, the FL is unsta- 
ble against an anisotropic superconducting (SC) state 
at least at T = K, provided that no disorder exists. 
When I*(iuji,q) is antiferromagnetic, the FL is unstable 
against a s ing let SC state. For example, it has been pro- 
posed [3^, HlJ that the condensation of cfy-wave Cooper 
pairs of the Gutzwiller quasiparticles due to the superex- 
change interaction is responsible for high-T c supercon- 
ductivity 32J, which occurs in the vicinity of the Mott 
M-I crossover. When I*{ioji,q) is ferromagnetic, the FL 
is unstable against a triplet SC state. 

The FL can also be unstable against a bond-order 
(BO) state, whose order parameter is { a \ a a j<y') and a 
flux state, which is a superposition of several BO states. 
Within Kondo-lattice theory, magnetic or SC states are 
more stable than BO and flux states. 

When U/W < 1, the conventional perturbation in 
terms of U is more useful than that in terms of I s (iu>i , q) . 
A similar argument to that for U/W > 1 applies in 
this case. When the nesting of the FS is sharp, a non- 
interacting electron gas is unstable against a spin density 
wave. When an interaction between electrons given by 
?7 2 Xs(»w;,q) is considered, the electron gas is unstable 
against an anisotropic SC state at least at T — K, 
provided that no disorder exists. 

No order parameter appears in one dimension. The 
FL for U /W > 1 becomes a Tomonaga-Luttinger liquid 
when I s (iu>i , q) is perturbatively treated, as the electron 
gas does when U is perturbatively treated. It is probable 
that Lieb and Wu's insulating state for the just half filling 
[33| can only be obtained by non-perturbative theory. 

In this analysis, the ground state is not a paramagnetic 
Mott insulator except for n = 1 and W/U = 0. How- 
ever, since the nature of actual M-I transitions cannot 
be explained within the Hubbard model, various effects 
should be considered in a multi-band or multi-orbital 
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model. Changes of lattice symmetries or jumps in lattice 
constants are often observed Q], which implies that the 
electron-phonon interaction should also be considered. 
A relevant electron-phonon interaction must arise from 
spin and orbital channels [34J. Cooperative J aim- Teller 
or orbital ordering must be responsible for the change of 
lattice symmetries. As well as the electron-phonon inter- 
action, the orbital-channel exchange interaction can play 
a role in the ordering. 

The FS sum rule holds for the Gutzwiller quasiparti- 
cles. The ordinary rule holds in the absence of an or- 
der parameter, and a modified rule holds even when the 
Brillouin zone is folded by an antiferromagnetic or or- 
bital order parameter. Then, Wilson's classification of 
crystalline solids into metals and insulators [IH can ap- 
ply to M-I transitions. Two types of M-I transitions are 
possible according to the band structures of the quasipar- 
ticles in the absence and presence of an order parameter: 
between a metal and an insulator and between a compen- 
sated metal and an insulator. The Kondo temperature 
Tjc can be different in metallic and insulating phases of 
a first-order M-I transition. In the metallic phase, Tk is 
higher than T and the quasiparticles are well defined. If 
Tk is lower than T in the insulating phase, the quasipar- 
ticles are not well defined. In such a case, the M-I transi- 
tion is a transition between a high-TV itinerant-electron 
phase and a Iow-Tr- local-moment phase. Change of lat- 
tice symmetries or jumps in lattice constants must play 
a crucial role in any first-order M-I transition. 

Since disorder, either small or large, must always exist, 
Anderson localization and other effects due to disorder 
[2^ can play a role in actual M-I transitions or crossovers. 

In conclusion, the ground state of the Hubbard model 
is not a paramagnetic Mott insulator except for n = 1 
and W/U = in two and higher dimensions. In order 
to explain actual M-I transitions, the electron-phonon 
interaction, multi-band or multi-orbital effects, and the 
effects of disorder should be considered beyond the Hub- 
bard model. Except for the Iow-Tr- local-moment phase 
of electrons, whether a crystalline solid is a metal or an 
insulator can be explained by an extended Wilson's clas- 
sification of the band structure of quasiparticles in the 
absence or presence of an order parameter. 
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